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1. Introduction 

Conservation laws are central to the analysis of physical field equations by providing 
conserved quantities, such as energy, momentum, and angular momentum. For a given 
field equation, local conservation laws are well-known to arise through multipliers pQ, 
analogous to integrating factors of ODEs [2 , with the product of the multiplier and 
the field equation being a total divergence expression. Such divergences correspond to 
a conserved current vector for solutions of the field equation whenever the multiplier is 
non-singular. If a field equation possesses a Lagrangian, Noether's theorem pQ shows 
that the multipliers for local conservation laws consist of symmetries of the field equation 
such that the action principle is invariant (to within a boundary term). Moreover, the 
variational relation between the Lagrangian and the field equation yields an explicit 
formula for the resulting conserved current vector. This characterization of multipliers 
for a Lagrangian field equation has a generalization to any field equation by means of 
adjoint-symmetries [Sill], whether or not a Lagrangian formulation exists. 

Recall, geometrically, symmetries are tangent vector fields on the solution space of 
a field equation and thus are determined as field variations satisfying the linearization of 
the field equation on its entire solution space. Adjoint-symmetries are defined to satisfy 
the adjoint equation of the symmetry determining equation on the solution space of 
a field equation |5 3 . (As such, unlike for symmetries, there is no obvious geometrical 
motion or invariance associated with adjoint-symmetries.) Through standard results 
in the calculus of variations pQ, it is known that the multipliers for local conservation 
laws are precisely adjoint-symmetries of the field equation subject to a certain adjoint 
invariance condition This allows a system of determining equations for multipliers 

to be formulated in terms of the adjoint-symmetry determining equation augmented by 
extra determining equations |H E]. In addition, the resulting conservation laws are 
yielded by means of a homotopy integral expression PH EH E] involving just the field 
equation and the multiplier, which is derived from the adjoint invariance condition, 
analogously to the line integral formula for first integrals of ODEs j2]. 

The purpose of this paper is to show that in the physically interesting situation 
where a field equation possesses a scaling symmetry, then the adjoint invariance 
condition and homotopy integral formula can be completely by-passed for obtaining 
conservation laws. In particular, a simple algebraic formula that directly generates 
conservation laws in terms of adjoint-symmetries for any such field equation is presented. 
Most important, when applied to a multiplier, the formula recovers the corresponding 
conservation law determined by the multiplier, to within a proportionality factor. This 
factor turns out to be the scaling weight of the conserved quantity defined from the 
conservation law. Consequently, all conserved quantities with non-zero scaling weight 
are obtainable from this formula [Hj. 

In Sec. the conservation law formula is derived. Examples and applications of 
this formula are presented in Sec. El As new results, first, a recursion formula is obtained 
for the local higher-order conservation laws of the sine-Gordon equation and a vector 
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generalization of the Korteweg-de Vries equation; second, a simple proof is given for 
closing a gap in the classification of local conservation laws of the Yang-Mills equations 
and Einstein gravity equations. Some concluding remarks are made in Sec. EJ 

2. Conservation Law Formula 

Consider a general system of field equations 

T A {x,u,du,d 2 u,...) = (2.1) 

for field variables u a (x) depending on a total of n > 2 time and space variables x a , with 
d k u denoting partial derivatives u a ^ ai ... ak = d k u a {x) / dx ai ■ ■ -dx ak , up to some finite 
differential order. (The coordinate indices a, [3, 7 run to n— 1; the field index a runs 
1 to N; the equation index A runs 1 to m. Summation is assumed over any repeated 
indices. This formalism allows the number of components of the fields u a and equations 
T A to be different.) For simplicity of presentation, the differential order of system (|2.1j) 
will be restricted to k < 2. 

Symmetries of the field equations (j2.1|) are the solutions 5u a = r) a (x, u, du, d 2 u, . . .) 
of the linearized equations 

C r ( V ) A = V a T A a + (D aV a )T Aa a + (D a D pV a )T A f = (2.2) 

for all u a (x) satisfying system (|2.1j) . where D a is the total derivative with respect to x a , 
and where T A a , T A ", etc. denote partial derivatives dT A /du a , dT A / du a )Q , etc.. The 
adjoint of equation ()2.2|) is given by 

C* T (uj) a = u A T A a - D a (u A T Aa a ) + D a D p (u; A T A f) = (2.3) 

whose solutions uj a (x, u, du, d 2 u, . . .) for all u a (x) satisfying system (|2.1|) are the adjoint- 
symmetries of the field equations (j2.1|) . Note the operators £t and C\ are related by 
the identity 

cuACriv f - r, a C* r {u) a = D a $ a {u, 77; T) (2.4) 

with 

$>, m T) = r/^AT A ; - D p (u A T Aa f)) - (D, V a )u; A T A f . (2.5) 

Hence, this expression f!2.5|) yields a local conservation law D a Q a (u), r/) = for any pair 
u A ,r] a , on all solutions of the field equations (|2.1j) . 

Now suppose the field equations are invariant under a scaling of the variables 

x a -> X pia) x a , u a -> X q(a) u\ (2.6) 

with p( a ^ = const, = const. From the corresponding scaling symmetry, given by 

S sU a = r] s a (x, u, du) = q {a) u a - p {a) x a u\ a) (2.7) 

the expression & a (u>, rj s ) produces a conserved current in terms of any adjoint-symmetry 
u A . 
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Proposition 2.1: For scaling invariant field equations (|2.1j) . every adjoint- 
symmetry ()2.3|) generates a conserved current on all solutions of ()2.1|) by the formula 

K = (q^-p^x^^T*;: - D p (u A T A f)) 

+ ((p(0 _ q i-) )u a p + pW^^^f. (2.8) 

Consider, now, a multiplier Q A (x, -u, 9m, d 2 u, . . .) for a local conservation law 

Q A T A = D^ a Q (2.9) 

that is assumed to be homogeneous ^Qj under the scaling symmetry, so 

5 S Q A = C Q (vs) A = r (A) Q A -p {a) x a D a Q A (2.10) 

with scaling weight r^) = const. Let s^ A > = const be the scaling weight of the field 
equations, 

5 S T A = C r (Vs) A = s (A) T A -p {a) x a D a T A . (2.11) 

Due to the scaling homogeneity of the constant r( A ) + s( A ' is independent of the 
index A. Then, the following important relation holds between the conserved currents 
*q ^d $J. 

Theorem 2.2: In terms of the scaling weights of the field equations (j2.11j) and the 
multiplier (|2.1U|) . every local conservation law (|2.9|) satisfies the scaling relation 

*5*u;q*5, w Q =r (A) + s^ + J2p {a) (2-12) 

a 

for all u a (x) satisfying the field equations, where "~" denotes equality to within a trivial 
conserved current ^T] DpG a P for some local expression 0°^ = — 0^°. Moreover, wq 
is simply the scaling weight of the flux integral J 'ty ( ^n a d n ~ 1 x defined on any (n — 1)- 
dimensional hypersurface x a n a = const (with normal vector n a ), 

5 S J VqU^x = w Q j %n a d^x. (2.13) 

Definition 2.3: A conservation law ()2.9)1 will be called (non)critical with respect 
to the scaling (j2.6J) if the scaling weight (j2.13J) of the corresponding conserved quantity 
is (non)zero. 

Corollary 2.4: As all multipliers necessarily are given by adjoint-symmetries, the 
conservation law formula ()2.8)1 consequently generates all noncritical conservation laws 
of the field equations (|2.1j) . 

The proof of relation ()2.12|) starts from the identity ()2.4|) with rj a = r] s a . We 
substitute the condition £* (T) a = — Cy(u) a on u A (holding for all u a (x) without use of 
the field equations), which is necessary and sufficient pQ for an adjoint-symmetry to be 
a multiplier Q A = u A . Here, £* is the adjoint of the linearization operator C w defined 
analogously to Cy and £y. We next use the adjoint relation 

Vs a C:(T) a = T A £ u (7k) A - D a ^ a (T, Vs ;uj). (2.14) 
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Substituting the scaling relations (|2.1(Jj) and (|2.11|) . followed by integrating by parts, we 
obtain 

w Q Q A T A = D a {®l + * a (T, Vs ; Q)) (2.15) 

where, note, the last term in this divergence vanishes when T A = 0. Then conservation 
law equation ()2.9|) leads to the scaling relation ()2.12j) . 

3. Examples and Applications 

3.1. Soliton equations 

For applications of the main conservation law formulas ()2.8|) and (j2.12|) . consider, firstly, 
soliton field equations in 1+1 dimensions. 

Korteweg-de Vries equation. The KdV equation in physical form for scalar field 
u{t, x) is given by 

T(u,u t ,u x , u xxx ) = u t + uu x + u xxx = 0. (3.1) 

This field equation is invariant under the scaling t — > A 3 t, x — > Xx, u — > \~ 2 u. The 
corresponding scaling symmetry 

5u = 1] S = —1u — 3tu t — xu x (3-2) 
is a solution of the linearized field equation 

£ r (v) = D t V + u x n + uD x n + D 3 x n = (3.3) 
for all u(t,x) satisfying the KdV equation. The adjoint of equation ()3.3j) is given by 

C* t (uj) = —D t uj - uD x uj - D 3 x u = (3.4) 

whose solutions uj for all u(t,x) satisfying (|3.1|) are the adjoint-symmetries of the KdV 
equation. Here, we see ^ Cr, reflecting the fact that the KdV equation (|3.1j) lacks 
a local Lagrangian formulation in terms of u(t,x). Now, for any adjoint-symmetry u, 
the conservation law formula ()2.8j) gives the conserved density 

$^ = -(3tu t + xu x + 2u)u. (3.5) 

From the obvious solution w = u of equation ()3.4|) . we consider the infinite sequence of 
KdV adjoint-symmetries ujq~) = (JZ*) k u, k = 0, 1, 2, . . ., generated by the operator 

n* = Dl + ~u+~D-\uD x ) (3.6) 

which is the adjoint of the well-known KdV recursion operator 

K = D 2 x + 2 -u+ l -u x D-\ (3.7) 

Note that, under the KdV scaling symmetry, u>^) — > X^ 1+k ^iU(k)- 

Each adjoint-symmetry uha is known to be a multiplier for a local conservation law 
of the form 

D t ^^(u,u x ,u xx , . . .) + D x ^ k} ( It, u x , u xx , . . . ) = (3.8) 
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on KdV solutions u(t,x), through lengthy calculations. For instance, originally the 
KdV conservation laws were derived one by one via the Miura transformation , from 
which the multipliers can be calculated. An alternative approach has involved extracting 
the conserved densities one by one through a residue method using a formal symmetry 
(pseudo-differential operator) for the KdV equation. More recently, in Ref. [Sj the 
conserved densities were obtained one at a time from a homotopy integral formula in 
terms of the adjoint-symmetries, after a verification of the adjoint invariance condition 
on each one. 

Here, by-passing such cumbersome steps, formula (|3.5|) yields the resulting 
conserved densities directly in terms of uq.), 

$ L - D x ( u xU(k)) - 2w^( fc) + 3t(uu x + u xxx )uJ(k) ~ -31l*(uj( k )) (3.9) 

which follows by means of properties of TZ*. This leads to a simple explicit recursion 
formula for all of the KdV local conservation laws 

Q 1 _|_ ok 

to within a trivial conserved density D X Q, as a result of the scaling formula (j2.12j) . 
Sine-Gordon equation. The sine-Gordon equation is given by 

T(u,u tx ) = u tx - sinu = (3.11) 

for scalar field u(t,x). This is a Lagrangian field equation with the scaling invariance 
t — > A _1 t, x — > Ax, u — > u. The corresponding scaling symmetry is 5u = r] s = tu t — xu x 
which is a solution of the linearized field equation 

C r {v) =D t D x r)-(cosu)r) = (3.12) 

for all u(t,x) satisfying the sine-Gordon equation. Here, as we have Cy = Cy, 
symmetries are the same as adjoint-symmetries, us = rj. Hence, for any symmetry 
rj, the conservation law formula ()2.8|) gives the conserved density 

®l = {xu x - tu t )D x j]. (3.13) 

We now consider the well-known infinite sequence of symmetries rf® = (lZ) h u x , 
k = 0,1,2, .. ., generated by the sine- Gordon recursion operator 

K = D 2 x + u x D- 1 (u x D x ), (3.14) 

starting from the translation symmetry rj = u x . Under the sine-Gordon scaling, note 
7/^ — ► \~ 2kr q^ k \ Each symmetry r/^ is a multiplier for a local conservation law of 
the form ()3.8j) on sine- Gordon solutions u(t,x). These conservation laws were found 
originally by an application of Noether's theorem ^3] and subsequently were derived 
by the same techniques used for the KdV equation [l], yielding the conserved densities 
one at a time through lengthy calculations. Here, similarly to the KdV case, formula 
()3.13|) directly leads instead to a simple explicit expression for all of the sine-Gordon 
local conservation laws, 

<^ ~ -D-\u x D xV {k) ) = -n nonloc (r] {k) )/u x = $J fc) (3.15) 
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and hence 

*;^Trk*w = rT2i(" <l+1) -^ w )^ < 3 - 16 > 

from the scaling formula (|2.12j) . Here 7£ non i oc stands for the nonlocal part of 1Z. If the 



relation — (D x & k ^) /u x = D x rf k ^ is substituted into expression (j3.16|) to get 



1 

then equation ()3.15j) yields an explicit conservation law recursion formula 



t/* +1 > = -u x *\ k) - D x (-D x *\ k) ) (3.17) 



+ 

with 

£ = u x 2 D x -^D x + ^ 2 + D~\C^ + ^m, 2 )/^) (3.19) 

representing a recursion operator on conserved densities. 

Modified Korteweg-de Vries vector equation. It is known that the recursion operator 
and symmetry hierarchy of the sine-Gordon equation are closely related to that of the 
modified Korteweg-de Vries equation u t + \u 2 u x + u xxx = 0. This scalar field equation 
has an interesting generalization jTH] 
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T(u, ut, u x , u xxx ) = u t + -u-u u x + u xxx = (3.20) 

for an N- dimensional vector field u(t,x), with any N > 1. The vector mKdV equation 
(13.20)1 is invariant under the scaling t — ► X 3 t, x —>■ Ax, w — > A _1 n. Its symmetries are 
the solutions of the linearized field equation 

3 3 

Cf(f)) = D t ff + 3u-f} u x + -u-u D x ff + D x ff = (3.21) 

for all mKdV solutions u(t,x), while its adjoint-symmetries are the solutions of the 
adjoint of equation ()3.21|1 

3^ _ _ _ ^ ^ ^ 3 ^ 

£*t(cJ) = —D t uj — -xu-u D X Q + 3u x -u> u — 3u-u x uo — D x uj = 0. (3.22) 

The vector mKdV equation ()3.2()j) admits the recursion operator 

U = D 2 x + u-u + u x D-\u- )-ujD-\u x A ), (3.23) 

where "A" denotes the antisymmetric outer product of two vectors and " j" denotes 
the interior product (i.e. contraction) of a vector against a tensor, namely c j(a A b) = 
(c • a)b — (c • b) a. This expression ()3.23|) is a manifestly SO(N) -invariant version of the 
vector mKdV recursion operator first derived in Ref. The adjoint of TZ is given by 
the similar operator 

TZ* = D 2 x + uj(uA )+u x ^D x \uA )+uD x \u-D x ) (3.24) 

(closely resembling the form of the sine-Gordon recursion operator in the scalar case 
N = 1, when the "A" terms vanish). We now consider the infinite sequence of mKdV 
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adjoint-symmetries generated by uji k \ = (7Z*) k u, k = 0,1,2,..., starting from the 
obvious solution of equation (j3.22j) . uj = u. By means of the scaling symmetry 

5u = ff s = — u — 3tu t — xu x , (3.25) 
the conservation law formula ()2.8|) yields the conserved density 

$ L = ~(3*Wt + + -D~^(u- D x u {k) ) = $' (fc) . (3.26) 

Note this can be expressed in terms of the recursions operator 1Z* by 

$ (fc) = ■ K iOB \oc& {k))/u-u x - (3.27) 

As urfy — > X^^ 1+2k ^cu^ under the mKdV scaling symmetry, we see from the scaling 
formula (j2.12j) that each adjoint-symmetry ujn^ is a multiplier for a local conservation 
law on mKdV solutions, given by 

D t *( fc )(«, Ux, u xx , ...) + D x ^ x {k) (u, u x , u xx , . . .) = (3.28) 

with 

*L ~ ~ $l v (3.29) 

However, in contrast to the scalar KdV and sine-Gordon cases, here expressions (j3.26J) 
and ()3.29|) do not lead in any immediate way to a recursion operator for the mKdV 
conservation laws (since D x u^ cannot be expressed directly in terms of $, fc x, w, and 
their derivatives). Nevertheless we have an explicit recursion formula 

*« = i^DZ\u-D x {{K*) k u)) (3.30) 

for generating all of the local conservation laws ()3.28|) (to within a trivial conserved 
density D X Q). 

Other soliton field equations, like the nonlinear Schrodinger equation, Tzetzeica 
equation, Harry-Dym equation, Boussinesq equation, and their variants ^7j, as well as 
more general multi-component scalar/vector field equations [HI], can be treated in a 
similar way to the preceeding examples. 

3.2. Fluid flow and wave propagation 

Secondly, field equations for fluid flow and nonlinear wave propagation in 2+1 and 3+1 
dimensions will be considered. 

Euler equations. The field equations for an incompressible inviscid fluid in two or 
three spatial dimensions are given by 

T(u, d t u, dga, P) = d t u l + ^djv! + p^&P = 0, T(d £ u) = 9 t V = (3.31) 

for fluid velocity u l (t,x) and pressure P(t,x), with constant density p. This system is 
invariant under the family of scalings t — > X p t, x l Xx l , u l — > A 1_p m*, P — > \ 2 ~ 2p P, 
for arbitrary p = const. The fluid symmetries are solutions (rf,r)) of the linearized field 
equations 

Cffort)* = D t rf + u j Djri + dju'if + p^D'n = 0, Cy {rff — Dirf — (3.32) 
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for all u(t,x), P(t,x) satisfying the Euler equations. The adjoint of equations (|3.32|) is 
given by 

C%(u,uj)i = —D t uji — u'DjijJi + diV?ujj — DiUj = 0, C\(ol) = —p~ 1 D l oOi = (3.33) 

whose solutions (u>i, uj) for all u(t, x), P(t, x) satisfying (|3.31j) are the adjoint-symmetries 
of the Euler equations. In addition to scaling symmetries, 

Su* = rjJ = (1 - P)u l - x j dju l - ptd t u\ (3.34) 
5P = r] s = (2 - 2p)P - x j djP - ptd t P, (3.35) 

the Euler equations are well-known to possess the Galilean group fHJ of symmetries, 
comprising time translations 

Su* = rf = dtv! = -^djv! - p^cfP, 5P = rj = d t P, (3.36) 
Galilean boosts, with velocity v % = const, 

8vf =if = titdju' -v\ 6P = r] = tv j djP, (3.37) 
and space translations and rotations 

Sv* = rf = C e u* = gdjU* - (vPdjg), 5P = ri = C^P = tfdjP, (3.38) 

where Ci x ) is a Killing vector of the Euclidean space in which the fluid flow takes place, 
i.e. =0. (In particular, £ l = a 1 — const yields translations, and £ l = Wxp 

b l i = —V 1 = const, yields rotations.) Corresponding adjoint-symmetries are given by 
the relations 

Ui = dif/d(tfu j ), (3.39) 
oo = J ooM - dttoidx 1 + Pdrf/d{d i P), (3.40) 

yielding 

oo, = 5 ijU j , oo = U tJ u % v? + p- l P, (3.41) 

-w w = ^V', (3.42) 

Ui = Sytv*, oo = S^itvP - x j ). (3.43) 
Now, typically, local conservation laws 

D t ¥(x, u, P) + D { ¥(x, u, P) = (3.44) 

on solutions of the Euler equations are derived through consideration of Newton's laws 
applied to fluid elements jTHj or by Noether's theorem in a Hamiltonian formulation 
Pl^]. In contrast, the conservation law formula ()2.8|) in terms of any adjoint-symmetry 
(oOi,oo) directly yields a conserved density 

= ((! - pW - rfdjrt - Ptdtu^ooi ~ ((ptu j - x j )d j u i + (1 - pypi (3.45) 

to within a trivial conserved density Z)jG\ Here, this formula easily leads to momentum 
and angular momentum 
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from the Killing vector adjoint-symmetries ([3.42 1) . and Galilean momentum 

= <J y toV ~ & 9 /w 9 (3.47) 
from the boost adjoint-symmetry ()3.43jl . as well as energy 

m l = -5-M ~ Ql/wv (3.48) 

from the fluid velocity adjoint-symmetry 1)3. 41 j) . to within proportionality factors. (A 
useful identity in these calculations is = D^x^u 1 ) on fluid solutions.) From the scaling 
formula ([2.12)1 . in three dimensions, it follows that = 4 — p when £ is a translation, 
wg — 5 — p when £ is a rotation, wn = 5 — 2p and = 4, representing the scaling 

weights of, respectively, the integrals for momentum and angular momentum J £-u d 3 x, 
energy f ||w| 2 d 3 x, and Galilean momentum J tv-u d 3 x. Note, for the dilation scaling 
p — 1, all the scaling weights are positive and hence the conservation laws (|3.46jl . (j3.47)l . 
(I3.48p are noncritical. These weights decrease by 1 in two dimensions, leading to the 
same conclusions. 

The Euler equations also are known to possess a vorticity conservation law 
|2Uj . which is unrelated to symmetries in contrast with the energy and momentum 
conservation laws ^Hj- Here, a derivation will be given by formula 1)3.45)1 directly in 
terms of fluid adjoint-symmetries. 

In three dimensions, vorticity is the curl of the fluid velocity 

= £ jk &u k (3.49) 

satisfying the vorticity equations [THj 

d£T + v?dff - Q'dj-u* = 0, 80 = 0, (3.50) 

where e t j k is the cross-product operator (i.e. e^ k = S h e l j k = is the totally 

antisymmetric symbol). We observe these equations have precisely the form of the 
adjoint-symmetry equations ()3.33j) and hence 

Ui = 5 {j n j , u = (3.51) 

yields a corresponding fluid adjoint-symmetry. Then the conserved density formula 
flUSt leads to 

^ = ^ (3-52) 



where, from scaling formula ()2.12jl . = 4 — 2p is the scaling weight of the vorticity 
integral J u- (d X u)d 3 x. Physically, this conserved quantity describes the total helicity 
(degree of knottedness) of vortex filaments. Note its scaling weight is noncritical 
provided p ^ 2. 

The situation in two dimensions is slightly different. The role of fluid vorticity is 
played by the scalar curl 

Q = e jk d j u k (3.53) 
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where e- k = £r. fe , is the antisymmetric symbol. This scalar vorticity satisfies the 
conservation equation 

d t Q + u j djVL = 0. (3.54) 

By taking a curl, we immediately see that 

Ui = e^&fL, uj = -tt 2 /2 (3.55) 

satisfy the fluid adjoint-symmetry equations (j3.33j) . The conserved density formula 
()3.45|) now yields 

*n = \^ ]k d ] u k ) 2 ~ tfjwn (3.56) 

with the proportionality factor u>q = 2 — 2p given by the scaling weight of the conserved 
vorticity integral f(d x u) 2 d 2 x (where d x u denotes the scalar curl e^ k d^u k ). Note 
that this vorticity quantity is noncritical if p ^ 1, i.e. other than for a dilation scaling. 
More generally, any function f{Q) in two dimensions is also a conserved density due 
to conservation ()3.54j) of the vorticity Q = d x u, but the resulting vorticity integral 
J f(d x u)d 2 x will have a well-defined scaling weight only if p — 0, namely for a spatial 
dilation scaling. In this case the conserved density 

** = f(e jk d>u k ) ~ & f /w f (3.57) 

again arises directly from formula (|3.45|) . through the adjoint-symmetry 

^ = e/n/", u = f-Qf. (3.58) 

Here the proportionality factor is simply Wf = 2, due to the scaling invariance of f(£l), 
and consequently the vorticity integral is noncritical. The same conclusion holds for the 
vorticity integral in three dimensions if a spatial dilation scaling, p = 0, is considered. 

The Navier-Stokes equations and polytropic gas dynamics equations can be treated 
analogously. 

Nonlinear wave equation. The scalar field equation for u(t,x) 

T (u, d t u, d s u) = d 2 u - &diU ± u u = -g aP d a d p u ± u v = (3.59) 

describes a nonlinear wave with interaction strength depending on a positive integer 
v > 1 (where g a p is the Minkowski metric tensor and x a = (t, x) are Minkowski 
spacetime coordinates). This is a Lagrangian wave equation invariant under the scaling 
x a — > Xx a , u — > X q u, for q = 2/(1 — v) ^ 0. The corresponding scaling symmetry 

5u = i] s = qu — x a d a u (3.60) 

is a solution of the linearized field equation 

C T (n) = -g al3 D a Dpu ± vu v - x i] = (3.61) 

for all u(t,x) satisfying the wave equation ()3.59|) . Here, symmetries are the same as 
adjoint-symmetries, since = Cy. The additional spacetime symmetries (Poincare 
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group) of the wave equation (J3.59)) for arbitrary v > are given by translations, 
rotations, and boosts, 

8u — Tj — C^u = £ a d a ii (3.62) 

where £ a (x) is a Killing vector of Minkowski space, d^ a ^ = 0. For certain interaction 
powers, v = 5 in 2+1 dimensions and v = 3 in 3+1 dimensions, the wave equation 
()3.59j) also admits [21] inversion symmetries 

1 

6* 



5u = rj = C^u + -div£ u (3.63) 



m a = -T a p{u, du)f - —^(div^ ud a u - -(d a div£)u 2 ) (3.65) 



T a p(u, du) = d a ud p u - -5 a ^ud lU ± —^u^ 1 ). (3.66) 



(where div£ = d a ^ a ) associated with conformal Killing vectors £° = c^x^x a — ^"x^x 13 , 
= const, satisfying d^ a ^ = Vtg a ^ for a conformal factor Vt(x). Through Noether's 
theorem, all spacetime symmetries ()3.62|) and ()3.63j) are known to be multipliers [21] 
for local conservation laws on solutions of the wave equation ()3.59|) . 

D a ^(t, x, u, d t u, d s u) = (3.64) 

with 

— - (div£ ud a u- -[ 
v + V s 2 

given in terms of the conserved stress-energy tensor 

1 2 

2 py 7 v + l 

The conservation d a T a /3 (u,du) = of this tensor ()3.66|) on solutions of the wave 
equation provides a well-known alternative derivation of the conservation laws associated 
with spacetime Killing vectors, 

#| = -T a p(u, du)£? for = 0. (3.67) 

Here, the resulting conserved densities instead will be obtained from the conservation 
law formula ()2.8|) directly in terms of the corresponding symmetries ()3.62|) and ()3.63|) . 
This yields 

$° = (( q - l)d a u - x p d a dfsu)(Cd 7 u + ^ydive u) 

- {qu - x^dpujid^^d^u + Cd a d^u + ^dw£d a u + — ^(<9 Q divO«) 

~ -w^l (3.68) 

to within a trivial conserved density (DpQ a/3 , with Q alS = —Q^ a ). The proportionality 
factor is, by the scaling formula (|2.12j) . the scaling weight of the flux integrals 
J ^r^dX, on a t = const spatial hypersurface S with normal vector n a = d a t. For 
translations C, a = a a = const, we have 

iy { = 4/(l-i;) (3.69) 

in 2+1 dimensions, while in 3+1 dimensions, 

w ( = {b-v)/(l-v); (3.70) 
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the weight increases by 1 for rotations and boosts £ a = b a @Xa, b alS = — b^ a = const, 
and increases by 1 again for inversions ()3.63|) so thus = 1 in the case of proper 
conformal Killing vectors. Hence, a critical case = only occurs for translations in 
3+1 dimensions when v — 5, and for rotations and boosts in 3+1 dimensions when v = 3 
as well as in 2+1 dimensions when v = 5, corresponding to scaling invariance of the 
energy-momentum integral J — T a ^{u, 9w)a^n Q ,dS and of the angular-boost momentum 
integral J —T a ^(u, dujb^x^n^T, in these cases. Consequently, $^ ~ is trivial only 
for these critical interaction powers [22] arid Killing vectors. In this situation, all local 
conservation laws are produced nevertheless from the more general formula (j2.5|) directly 
in terms of pairs of Killing vector symmetries (|3.fi2j) , 

y a ( m , m )^-T%(u,du)f (3.71) 

where £ a = f^,^]" ^ s the commutator of the Killing vectors. The same result holds 
even for the non-critical cases where ^ 0. 

Other nonlinear wave equations, such as sigma models and wavemap equations j2H|, 
can be treated in the same way. 



3.3. Gauge theories 

Finally, Yang- Mills fields and gravitational fields in 3+1 dimensions will be considered. 

Yang-Mills theory. The Yang-Mills field on Minkowski space (IR 4 ,^^) is a vector 
potential A a a (x) that takes values in an internal Lie algebra Q = (R , c°f, c ). Associated 
with A a a is the Yang-Mills covariant derivative 

V£ = d a + c\ c A\ (3.72) 

and the Yang-Mills field strength tensor 

Kfi = + \^A b a A c (3.73) 

where c a b c denotes the structure constants of the Lie algebra Q. The Yang-Mills equation 
with gauge group based on Q is then given by 

r;(A, dA, d 2 A) = g^F^ = (3.74) 

which is invariant under the scaling x a — > Xx a , A a a — > A _1 A^. Whenever the gauge 
group is semisimple, the Yang-Mills equation ()3.74j) arises from a Lagrangian (see e.g. 
Ref. |211), an d i n this situation both the symmetries and adjoint-symmetries of this field 
equation (J3.74j) are given by solutions 5A a a = r\ a a of the linearized Yang-Mills equation 

£r(v); = ff^Vy^ + c\ cV b a F^) = (3.75) 

for all Yang-Mills solutions A a a (x). Note, here, = D a + c a ^ c A h a acts as a total 
derivative operator. The well-known local symmetries of the Yang-Mills equation ()3.74j) 
are comprised by gauge symmetries 

<K = V a X° (3.76) 
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involving any Lie-algebra valued scalar function x a { x i A 9 A, . . .), and spacetime 
symmetries 

5A a = = _ V^AJ) (3.77) 

where q(x) is any conformal Killing vector on Minkowski space, 

the case of a dilation Killing vector, $ — x@, the spacetime symmetry ()3.77|) reduces to 
a sum of the Yang-Mills scaling symmetry and a gauge symmetry, 

= < = ~x%A a a -A a a = -2x?F% - V A a {x p A a p ). (3.78) 

A recent classification analysis [22] has proved that these are in fact the only nontrivial 
local symmetries admitted by the Yang-Mills equation if the Lie algebra Q is real and 
simple. However, for a simple Lie algebra Q with a complex structure, the same analysis 
found that the Yang-Mills equation also admits complexified spacetime symmetries 

&< = 2j\fF b pa (3.79) 

where j a b is the complex structure map on Q, satisfying the properties 

jVc = -€ 3 a b c b cd = c\ d f c , k a[b f c] = (3.80) 

with k ab = c c a dC d bc being the Cartan-Killing metric on Q. Through Noether's 
theorem, these symmetries ()3.76|) to (|3.79j) are multipliers for local conservation laws 
D a ^/ a (x, A, dA) = consisting of, respectively, 

= 9 av 9^k ab Fl^ X b ~ (3.81) 

related to the Bianchi identity on F^, and 

q a = T a p(F)f, ^ a = f a p {F)^, (3.82) 

given by the conserved Yang-Mills stress-energy tensor 

T aP {F) = g u °k ab (FlFl - \g a ^Fl v F b la ) (3.83) 
and its complexification 

T aP {F) = g^J ab (K^ - \s^F% F \°)- ( 3 - 84 ) 

These conservation laws ()3.82j) yield (complexified) energy-momentum, angular and 
boost momentum for translation, rotation and boost Killing vectors = 0, and 

additional quantities for dilation and inversion Killing vectors cK Q £^ = |div£g Q/3 ^ 0. 

However, the previous results do not fully settle the classification of local 
conservation laws of the Yang-Mills equation ()3.74|) . since it leaves open the question 
of whether any trivial symmetries could yield nontrivial conservation laws by Noether's 
theorem. For Lagrangian field equations whose principal part (i.e. highest derivative 
terms) is nondegenerate, it is known that there is a one-to-one correspondence between 
nontrivial variational symmetries and nontrivial conservation laws jU Ej. But this 
correspondence is not automatic for a field equation with gauge symmetries, due to the 
resulting degeneracy of the field equation's principal part, in contrast to the previous 
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examples in this section. Here, through an application of formula ([2.8)1 . the gap in the 
classification of Yang-Mills conservation laws will be addressed. 
We consider local symmetries 

5A a a = v a a (x,A,dA...) (3.85) 

assumed to be homogeneous with respect to the Yang-Mills scaling ()3.78|) . 

S B V a a = r vl - ^d P 7] a a (3.86) 

with scaling weight r = const. As noted in Sec. El there is no loss of generality in such a 
homogeneity restriction. Now, formula ()2.8)1 yields a conserved current generated from 
any such symmetry, 

$ a (?7, S S A) ~ 2k a \x^Fl p g p[v V Aa] r 1 b v - rj a u (x^F bav + F bau )), (3.87) 

with the current being linear and homogeneous in 77^ and V^r/^. If 77° = Q a a is a 
multiplier for a local conservation law of the Yang-Mills equation (|3.74jl . 

K b 9 a "Q> Av Kp = D ^Q ( 3 - 88 ) 
where the current 

^q{x,A,8A,...) (3.89) 

can be assumed homogeneous under the Yang-Mills scaling ()3.78|1 . then the scaling 
formula ()2.12|) gives the relation 

® a {r),5 s A)~(r + l)% (3.90) 

to within a trivial conserved current. Hence, for a variational symmetry that is trivial, 
so 77° = for solutions of ()3.74|) . we see that, if r 7^ — 1, 

^^^$ a (0,M)=0 (3.91) 

is a trivial noncritical current on Yang-Mills solutions. Note the factor r + 1 here is 
precisely the scaling weight of the flux integral of the current This establishes the 
following classification result. 

Proposition 3.1: No nontrivial conservation laws ()3.88j) that are noncritical — i.e. 
whose associated conserved quantity f ^/ a d a tdT,, on a spatial hypersurface £ given by 
t = const, has non-zero scaling weight — may arise from variational trivial symmetries 
(j3~%5) and (j3~£fij) of the Yang-Mills equation (jHZZSj) - 

General Relativity. The vacuum gravitational field equation on a 4-dimensional 
spacetime manifold is given by 

T a ^g,dg,d 2 g) = G a ^ = (3.92) 

where G a/3 is the Einstein tensor (i.e. trace- reversed Ricci tensor R a/3 ) for the spacetime 
metric g a p{x). Its linearized field equation is given by 

C r {r]) a(3 = -i(V^V^ Q/3 - 2V^V 9{ V )/J + 9 aP ^lv^) = (3.93) 
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whose solutions for all g a g{x) satisfying the Einstein equation (|3.92j) are the symmetries 
5g al3 = i] af3 (x, g, dg, . . .) of the gravitational field, where is the covariant total 
derivative operator associated with the metric, V 9 ^g a/3 = 0, and "bar" denotes trace- 
reversal on a given tensor. Here, since the gravitational field equation (|3.92|) comes from 
a Lagrangian adjoint-symmetries are the same as symmetries. It is known that the 
only admitted nontrivial local symmetries [2Z1 consist of a constant conformal scaling 

S9 a p = 9 a p (3-94) 
and diffeomorphism gauge symmetries 

89ap = C C9 a p = ^\ a C p) (3-95) 

for any local vector field (^(x, g,<9g, . . .). If we consider local conservation laws of the 
gravitational field equation ()3.92|) . 

V 9 a * a (x,g,dg,...)=0, (3.96) 

then through Noether's theorem the constant conformal scaling 1)3.94)1 is not a multiplier, 
while the diffeomorphisms ()3.95j) yield a trivial conservation law = (^G^ = on 
solutions of ([UnSD - 

Ref. [2Zj asserts that in fact there exist no nontrivial local conservation laws ()3.96)) . 
but does not provide a full statement of the proof. Here a complete classification proof 
will be given for diffeomorphism-covariant conservation laws, by use of a covariant 
version of the conservation law formula f)2.8j) . A suitable dilation scaling symmetry 
is provided by the diffeomorphism symmetry (J3.95)) specialized to a homothetic vector 
field ( a = £ a , namely 

5 a g af} = Csg, V 9 J a = const ^ 0. (3.97) 

Now, a local conservation law ()3.96|) is diffeomorphism-covariant whenever it is of the 
form 

V 9 Jf a (g,R,V g R,...) = (3.98) 

on solutions of the Einstein equation ()3.92j) . with satisfying the natural 
transformation property 5^ a = C^ a under all diffeomorphism symmetries ()3.95|) . 
Correspondingly, without loss of generality we consider diffeomorphism-covariant local 
symmetries 

Sg aP = r ] a ^g,R,V 9 R,...), (3.99) 
which are necessarily homogeneous with respect to the dilation scaling 

<W^ = £^ a/3 - (3-100) 
Then for any such symmetry ()3.99|) . formula ()2.8|) yields a conserved current 

$ a ( v , C t g) = -(CR^ a ~ V* V^)^ + V^V^ + V^V^rT (3.101) 
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whose dependence on r] al3 and ^7 9 f j] al3 is linear, homogeneous. It now follows from the 
scaling formula (|2.12j) in covariant form that if rj = Q a/3 is a multiplier for a local 
conserved current, Q a gG a ^ = V^q, then to within a trivial conserved current, 

$ a (r),£sg)~w*«, w = V^V0. (3.102) 
Hence, for a trivial multiplier that is diffeomorphism-covariant ()3.99|) . we have 

tyq ~ uT^O, C^g) = (3.103) 

on solutions of the Einstein equation 1)3.92)1 . since 7] al3 = 0. Therefore, the following 
result holds. 

Proposition 3.2: No nontrivial diffeomorphism-covariant conservation laws (J3.98)) 
may arise from variational trivial symmetries ()3.99|) of the gravitational field equation 

Solutions of the Einstein equation ()3.92|) with a Killing vector C^g a/3 = 2^7^^ = 
are known to possess an unexpected extra symmetry [2H] 

^ aP = vg^-1i {aXpY 8t a = -xJi% (3-104) 
where v is the scalar twist and Xb is the dual of £ , defined by the equations 

V 9 ^ = e^f, 2Vf oX/3] = £^VT. (3-105) 

Note this symmetry ()3.104|) is nonlocal due to its dependence on v, Xg (in terms of 
g a p,£ a )- Here, a corresponding nonlocal conservation law will be derived through the 
conservation law formula ()2.8j) in covariant form using the constant conformal scaling 

® a (v,g) = 0„„v«V" - vy*. (3-106) 

This yields, from ()3.104|) . 

S"^) = + 2V^(e ( V } ) * 2^V^x a) - (3.107) 

Moreover, on solutions of the Einstein equation ()3.92|) with two commuting Killing 
vectors an d C2 > the symmetries 1)3.104)1 are known to generate an infinite-dimensional 
algebra of nonlocal symmetries (the Geroch group |29j). Then the covariant formula 
(13.106)) leads to a corresponding infinite sequence of conservation laws, related to the 
complete integrability of the system G af3 = 0, C^g a g = C^ 2 g a/3 = 0. 

A similar classification treatment of local and nonlocal conservation laws of the 
self-dual Einstein equation and self-dual Yang-Mills equation will be given elsewhere. 

4. Conclusion 

The conservation law expression derived in Proposition 2.1 is a generalization of a 
formula mentioned in Ref. PQ in the case of linear PDEs and extends a similar formula 
considered for self-adjoint PDEs in Ref. A variant of this expression has been 
central to a recently obtained classification of local conserved currents for linear massless 



Conservation laws of scaling-invariant field equations 



18 



spinorial field equations of spin s > [SOI EI]- Numerous applications to nonlinear ODEs 
are presented in Ref. [2|. 

Apart from its main use in generating all noncritical local conservation laws in an 
algebraic manner — by-passing the standard homotopy integral formula and (adjoint-) 
invariance conditions — in terms of local (adjoint-) symmetries admitted by any given 
PDEs, the conservation law expression is able to produce nonlocal conservation laws 
from any admitted nonlocal (adjoint-) symmetries. Such examples and applications will 
be considered in a forthcoming paper [32] . 
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